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Abstract 

The visualization of any graph plays important role in various aspects, such as graph drawing software. 
Complex systems (like large databases or networks) that have a graph structure should be properly 
visualized in order to avoid obfuscation. One way to provide an aesthetic improvement to a graph 
visualization is to apply a force-directed drawing algorithm to it. This method, that emerged in the 
60's views graphs as spring systems that exert forces (repulsive or attractive) to the nodes. The first 
force-directed drawing algorithm of Tutte [1] was aiming at providing a barycentric representation of 
the graph and it is suitable for getting a straight-line drawing of 3-connected planar graphs without 
crossings. Each vertex is gradually placed at the berycenter of its neighbour vertices. This method used 
no spring systems, though it can be easily understood that springs can be applied to the system getting 
the same result as that of the barycentric method. Unfortunately, Tutte's method usually gives bad 
angular resolutions. 

Many other force-directed methods have appeared since then. Eades' method [5] creates 2D layouts 
using spring systems and it is suitable for graphs with less than 30 nodes. Attractive forces are applied 
to adjacent vertices only and repulsive to all pairs of vertices (so we can view the edges as springs with 
zero rest length and the nodes as electrically charged identical particles). Of course, the exact type of the 
forces exerted is a very important aesthetic criterion. Eades sets logarithmic attractive forces and inverse 
square root repulsive forces. Fruchterman and Reingold [5] apply the notion of simulated annealing, in 
order to provide an even better 2D layout of less cumulative spring energy, compared to that of Eades. 
The main disadvantage is, again, that the algorithm is limited to graphs having less than 40 vertices. 
Hadany and Harel |4] provide a force-directed layout for much larger graphs of over 1000 vertices. In 
order to obtain that, they follow a multi-scale approach where they initially consider an abstraction of 
the graph, which is iteratively augmented by the graph details. Gajer, Goodrich and Kobourov [5] use 
a coarsening strategy which avoids quadratic space and time complexities. They use a vertex filtration 
which enables the algorithm to restrict the number of vertices relocated. They, also, introduce the idea 
of realizing the graph in high-dimensional Euclidean space. An excellent survey on various force-directed 
methods can be found in [6]. 

A Lombardi drawing of a graph is a drawing where the edges are drawn as circular arcs (straight edges 
are considered degenerate circular ones) with perfect angular resolution. This means, that consecutive 
edges around a vertex are equally spaced around it. In other words, each ang i<S between the tangents of 
two consecutive edges is equal to 27r/(i where d is the degree of that specific vertex. This kind of drawing 
took its name from the American artist Mark Lombardi, who mainly used this technique in his paintings. 
The requirement of using circular edges in graphs when we want to provide perfect angular resolution is 
necessary, since even cycle graphs cannot be drawn with straight edges when perfect angular resolution 
is needed. 

In this survey, we provide an algorithm that takes as input a random drawing of a graph and provides 
its Lombardi drawing. 



^ every angle in this survey will be expressed in rads 
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Figure 1. The spring forces on the central vertex due to the connection with the 4 circular 
edges-springs, form a cumulative force equal to the sum of the spring forces of 4 straight spring-edges 
(dotted lines), each one having the same length as that of its respective circular arc. In the specific 
example, however, we can easily see that the cumulative force on the central vertex is not zero, 

something which means that our algorithm will move the central vertex to a position where X]i=i 



1 Lombardi Graphs 

We will use the spring-force layout of a graph. A very important property of a single vertex which is 
connected to a number of d vertices by circular edges is that, if the edges are regarded as circular springs, 
the total spring force exterted on the central vertex will be equal to the one exerted on it if the circular 
springs are replaced by straight ones, each one having the same length as the respective circular arc and 
a direction onto the tangent of each respective circular edge (as in figure 1). The reason behind this, is 
that each force is exerted instantly on the direction of each tangent. In figure 1 there is an example of 
a node having degree equal to 4. Whether we take the 4 circular edges or the 4 straight ones (of equal 
corresponding length) under account, the instant result on the move of the central vertex will be the 
same. 

We will use this relation extensively while forming the angular resolution around a vertex. Our 
procedure actually uses, initially, an already known force-directed method and, afterwards, starts fixing 
the angular resolution around each vertex of the current drawing. The target is to obtain a drawing with 
the least cumulative spring energy, which will, also, have perfect angular resolution. 

The type of force-directed method initially applied, depends on the number of nodes of G. Let's 
assume that this is not great enough (at most 40 vertices), so we can apply the barycentric method of 
Tutte. However, there are methods, as predescribed, for much greater numbers of nodes. After applying 
the barycentric method, we order the nodes of the graph according to a descending order of degree. We 
start by picking the first vertex in the order and choosing its neighbour vertex that is farthest away in 
straight distance from i lIl(for example we pick vertex b in figure 2). We erase the current edge (circular 
or straight) between them and connect them with a straight edge (so in a case of a previous straight 
edge connection, we do not change anything). We continue in a clockwise order to the next adjacent 
edge (that is c in the example), which we erase and connect the 2 corresponding vertices (the central 
and its neighbour vertex) with a circular edge which has angle 2n/d with the previously formed edge 



■^the reason behind this is that the minimum total energy that we want to achieve forces us to connect initially the 2 
farthest points in plane with a straight edge; any other edge gives greater length and this means greater energy 
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Figure 2. The central vertex a has degree 4, so starting with neighbour vertex b and continuing in 
clockwise order, we arrange each new circular edge having angle equal to 7r/2 with the exactly previous 
one, until we reach again b where we stop. Bold (straight) edges refer to the previous edge visualization. 



(in anticlockwise order and we obviously refer to angles between tangents in the position of the central 
vertex). This is totally valid since any two points in plane can be connected with one circular arc, if we 
specify a tangent passing through one of them. We continue doing this procedure, creating equal angles 
between tangents, until we reach the very first edge, where we stop. 

It is obvious, now, that the Lombardi property has been applied on the central examined vertex. 
However, the zero cumulative spring force has been lost because the spring forces have changed, in both 
the direction and the length of the springs. So, we must fix that, by finding the expression of a spring 
force on a vertex when the spring is circular. Here is where the first observation is used. We must bring 
the central vertex to the position of zero cumulative spring force, leaving stable its neighbour vertices. 
So, our target is to minimize the energy of the spring system of this examined central vertex which is 
connected to d vertices under the invariance of the Lombardi property, which means that we are not 
allowed to alter the angles between the tangents of the edges and the x-axis, formed in the plane. In 
other words we seek to minimize the sum of the circular springs' lengths which equals where ti is 

the arc length of the i-th adjacent edge (starting from the one we beginned with in the previous step and 
continuing in clockwise ordeiH). When we accomplice that and move the position of the central vertex 
to a position of mimimum total spring energy, we continue with the rest of the vertices in the descending 
order and then we repeat the whole algorithm from the beginning several times to gain convergence. The 
number of times may vary for each drawing, so we can just calculate the difference in the co-ordinates 
of every vertex between each 2 repetitions of the algorithm and stop where we have reached the desired 
accuracy, for each one of the vertices. 



''so in our example ti is the circular edge (straight edges are regarded as circular arcs of infinite radius) connecting a 
with vertex b 



Figure 3. For any given arc length t, the corresponding chord length obeys the relation 

X — t sin{9)/6 where > is the angle between x and the tangent of the arc in the intersection point 

between the chord and the arc. 

2 The Algorithm 

A basic trigonometric relation is that the length of a chord x is equal to 

X = t* sin{d)/9 (1) 

where t is the length of the arc that is defined between the 2 points of the chord and 7r/2>6'>0is the 
angle between the chord x and the tangent at any of the 2 intersection points of the defined cycle and 
the chord, as we can see in Figure 3. 

The Euclidean distance between any points {pi,p2) in plane (where pi ~ {xi,yi), p2 = (2:2,1/2)) is 
given by the relation 

Dipi,P2) = - 2:2)2 + (2/1 - 2/2)2 (2) 

Since we seek minimization of the sum of arc lengths, while keeping the Lombardi property, we want 
to minimize the quantity 

d 

'^Xi* di/sin{9i) (3) 
1=1 

where Xi refers to the distance between the central vertex and the i-th neighbour vertex (assuming a 
clockwise enumeration as indicated before from i = 1 to d) and 9i is the angle between the tangent of a 
circular edge connecting the central vertex with the i-th vertex neighbour of it and the straight distance 
between the central vertex and the currently examined neighbour of it (0 < ^ < 7r/2). In order to avoid 
obfuscation with the co-ordinates we will refer to the distance Xi as D from now on, associated with each 
specific central examined vertex. 

From (2), (3) we get: 

d 

X] * y(2?r^^a;^^FH~(yi~^^/s*'^(6'0 (4) 

4=1 

The Lombardi property of equally spaced tangents on arc edges must be preserved while moving the 

central vertex on plane. In other words, we move the central vertex trying to obtain X]i=i where Fi 
refers to the spring force exerted between the central vertex and the i-th neighbour vertex of it, assuming 
zero rest length on each spring. 

However, sin{9i) is equal to the distance between the i-th neighbour node's (stable) position and 
the tangent of the specific arc (see figure 4) divided by the new Euclidean distance between the 2 ver- 
tices (new position of central vertex and the stable position of its neighbour vertex) which is equal to 



□Id position 



previous (straight) edge 



Figure 4. The central vertex is moved around trying to minimize the total cimulative spring force 
exterted on it and preserving the Lombardi property. Here we can see the central vertex moving to a 
new position along with one of its neighbours which stays stable. Of course the previous edge may have 
also been circular and not straight as here, which means that we should have used its tangent. 





(»,].yo) 







Figure 5. The distance between point {xa, Ha) and to the straight line defined on the line segment 
connecting points {xo,yo) and (a;h,?/h). 



{xi — Xo)"^ + (j/i — J/o)^, where the co-ordinates {xo, yo) refer to the new position of the central vertex. 
Assume that the central node's initial position is {xQ,yQ). In order to avoid complicated equations we 
think of the tangent of the first formed edge as tantamount to the semiaxis Ox and since we preserve the 
Lombardi property while moving around the central vertex, all the tangents keep the same ordering in 
plane (e.g. each one moves in parallel to its previous arrangement). This fact helps us constructing the 
arrangement of the d tangents if we set the first one onto semiaxis Ox where we think of the new position 
of the central vertex on (0,0). Therefore, we obtain ?/ = for the first tangent{footnotewhich we call initial 
from now on, refering to each specific central examined vertex, y = tan(— (z — l)*27r/d)x = tan(—2*TT/d)x 
for the second one (by setting i — 2) and so on, until we get to the d-th tangent. 

The only thing remaining is to identify the distance between each tangent and the corresponding 
neighbour vertex. In general, the distance between a straight line (identified by 2 points) and a point in 
plane, is given by the following equation, where the co-ordinates refer to figure 5: 

^ ^ (ajj, - Xa){yo - y^) - {xp - Xa){yb - yg) 
y/{xb - Xf))'^ + {yb - yaY 

In our case and since we define each tangent on a semiaxis starting at point (0,0) wc surely know point 
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i^OyUo)- Point {xa^ya) is the actual neighbour vertex i correlated with the specific tangent and point 
{^b,yb) can be obtained by setting xi, = I and using the tangent form {yb/xb) = tan{—(i — 1) * 2-K/d) 
for the specific i, to obtain yb- The overall distance between the central node's new position and the i-th 
neighbour node's (stable) position is, then: 

^ ^ \-yi+Xi* tan{-{i - 1) * 27r/d)| 
^yi + tan^{{i - l)2n/d) 

Using (6), we get for the sin{6i): 

sin{e,) = ^ (7) 

V {xi - xqY + {yi ~ y^Y 

Replacing (7) into (4), we take a relation solely depending on the central node's and i-th neighbour's 
coordinates. The sum of it, according to all neighbour vertices, can be minimized, moving the central 
vertex to a position of minimum total energy: 



d D * arcsin( 



D 



D (8) 

i=l 

which is equal to 



d 

min 



xo,yo 

i=l 



y^iixt - xoY + iVi - yo f] * arcsin{—j= ^ ) (9) 

^J[xi - xnY + {yi - yo) 



We note that D is defined in (6) and that there is no problem with the range of values of the arcsine 
function since < 9 < Tr/2. 

After altering the position of the central vertex and moving it to a balancing position, we apply the 
exact same procedure on the next vertex of the descending order and so on until we reach the |y|-th 
vertex in the order. There, we have to repeat the whole procedure from the beginning since the positions 
of the vertices should have changed, as many times as it needs for convergence under a desired accuracy. 

Finally, we must note that each repetition of the algorithm on a different central vertex should be 
coming along a rotation of the Cartesian coordinates in order to place the straight distance between the 
old position of the central vertex and the farthest edge onto semiaxis Ox, as we have described before. 
So, each time we switch to another central vertex in the descending order we should place the first 
tangent onto the new Ox semiaxis by multiplying the coefficients of all tangents (the first point needed 
to describe them is always (0,0) regarding all of them which is not rotated obviously, so we only apply 
the transformation on each second point) by 

cos{(j)) —sin{(j}) 
sin{(j)) cos{(j)) 

which corresponds to the clockwise rotation matrix of the 2-dimensional Cartesian coordinates, where 
angle is the angle between the 2 semiaxes (intersected at (0,0) of the two initial tangents of any 2 
consecutive examined central vertices). 
The pseudocode follows: 



repeat 

Apply a barycentric method to G=(V,E) 

Arrange the vertices of the graph according to descending order of degree; 
Let u = |F| and Vi be the ordering; 
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//start fixing tangents 

/ /there are no tangents in the beginning; these are formed in the first iteration of the FOR proeedure 

/ /and their direetion stays unaUered unless we unmark them, 
for i = 1 to V do 

Find the neighbour vertex of Vi which is farthest away in straight distance from it; 

Delete current edge between them and connect them with a straight edge, marking it with a 

tangent intersecting the central vertex; 

Change (Rotate) Cartesian co-ordinates in order to place this (initial) tangent onto semiaxis Ox; 
Get the next edge of the central vertex in clockwise order; 
start applying the Lombardi property around the central vertex 
for j = 1 to d-1 do 

Delete current edge between Vi and j and connect them with a circular edge having angle 2j7r/c? 
with the initial tangent; 

Mark the tangent on this circular edge intersecting the central vertex; 

j ^ j + 1 
end for 

Move the central vertex to a position of zero cumulative force on it according to equation (8) by 
not altering the directions of marked tangents in plane; 

When the position of the central vertex gives zero cumulative force according to (8), set i ^ i + 1 
and unmark all tangents; 
end for 

until the difference on the co-ordinates of each vertex converges 

3 Conclusions 

We have sketched an algorithm for creating a Lombardi drawing of a graph, given a random drawing 
of the graph. However, the algorithm's convergence remains to be tested for various numbers of graph 
nodes, since force-directed methods, in general, experience different convergence attitudes, which are 
highly correlated with the number of nodes of the graphs they are applied on. Especially in our case, 
graphs with very high degrees on vertices should be checked for appropriate convergence. 
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